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We investigate torque and also conventionally defined spin-Hall currents in two-dimensional 
(2D) spin-orbit coupled systems of spin-1/2 particles within the linear response Kubo formal- 
ism. We obtain some interesting relations between the conventional and torque spin-Hall con- 
ductivities for the generic effective Hamiltonian Hq — e% + A{}<i)ax — B(k)ay, where j4(k) = 
rjfki + r){jkikj + rjf^ikikjki + ■ ■ ■ , B(k) = rjfki + rjfjkikj + r)fjikikjki + ■■■ , and »7's are the spe- 
cific system-dependent coefficients. Specifically, we find that in the intrinsic case the magnitude 
of torque spin-Hall conductivity (7jj(0) is always twice larger than the conventional spin-Hall con- 
ductivity (T^j(O), and the two conductivities have the opposite signs, i.e., o-Jj(O) = —2(7%y{0). This 
universal relation, therefore, suggests that in the intrinsic case, the total spin Hall conductivity 
o"|y(0) in the 2D systems is equal to conventional spin Hall conductivity in magnitude but has the 
opposite sign, namely, <Jxy{0) ~ f"ly(0) + cr°j(0) = —a^y{0). This universal relation also holds 
in the presence of an uniform in-plane magnetic field. We also find that if the 2D systems are 
rotationally invariant, there exists a hyper-angular momentum = (k x §f)^Sz + Lz which is 
conserved. Furthermore, the hyper-angular momentum current {^{Iz,Vx}) vanishes, and this leads 
to a hyper selection rule for the conventional spin-Hall current. In particular, in the 2D k- linear 
Rashba and wurtzite-type systems, 1^ = Sz + L^, and the up(down)-spin current would always 
be accompanied by the down(up)-orbital angular momentum current (OAM). In the 2D k-cubic 
Rashba, Iz = Ssz + Lz, and the hyper-selection rule is the same as in the k-linear Rashba system. 
In the 2D k-linear Dresselhaus system, on the other hand, Iz ~ —Sz + Lz, and the up(down)-spin 
current would always be followed by the up(down)-OAM current. 

PACS numbers: 71.70.Ej, 72. 25. Do, 73.63.Hs, 85.75.-d 



I. INTRODUCTION 



perature has been detected electrically. [Ill d) Ell 



Spin current generation is an important issue in the 
emerging spintronics.[l|, 0, Q Recent proposals of the in- 
trinsic spin Hall effect are therefore remarkable [1, 
In the spin Hall effect (SHE), a transverse spin current is 
generated in response to an electric field in a system with 
spin-orbit coupling, 0] This effect has been consid- 
ered to arise extrinsically, i.e., by impurity scattering f6l|. 
The scattering becomes spin-dependent in the presence 
of spin-orbit coupling, and this gives rise to the SHE. In 
the recent proposals, in contrast, the SHE could arise in- 
trinsically in hole-doped (p-type) bulk semiconductors Q 
and also in electron-doped (n-type) semiconductor het- 
erostructures due to intrinsic spin-orbit coupling in 
the band structure. This intrinsic SHE would thus pro- 
vide a mechanism to generate electric driven spin current 
without applied magnetic fields in semiconductors, which 
can be more readily integrated with well-developed semi- 
conductor electronics. Recently, the spin accumulation at 
the edges of semiconductor samples which is believed to 
be due to the SHE, has been measured opticallvfsi. lol. ITdj . 
Further, large SHE in metallic systems even at room tem- 
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Many theoretical papers have been written addressing 
various issues about the intrinsic SHE. In Ref. the SHE 
in the p-type GaAs semiconductor was explained as aris- 
ing from the k-space Berry curvature in response to the 
applied electric field. This intrinsic SHE would lead to 
the possibility that the spin-orbit coupling can be used 
to manipulate spin chirality in semiconductors without 
dissipation. In Ref. it was shown that the SHE in 
the p-type GaAs semiconductor is robust against the dis- 
order based on the parity invariance of the spherical Lut- 
tinger Hamiltonian. The Berry-phase-induced SHE was 
also generalized to the case of spinning particles jT5|. In 
p^ . an orbital-angular- momentum (OAM) Hall current 
is predicted to exist in response to an electric field and is 
found to cancel exactly the spin Hall current in the SHE. 
In [17|, however, ah inito relativistic band structure cal- 
culations show that the OAM Hall conductivity in p-type 
semiconductors is one order of magnitude smaller than 
the spin Hall conductivity, indicating no cancellation be- 
tween the spin and OAM Hall effects in bulk semiconduc- 
tors. The spin Hall conductivity in the two-dimensional 
(2D) k-linear Rashba system has been shown to be sup- 
pressed by weak non- magnetic disorder [T8| . However, 
the spin-Hall conductivity calculated with the considera- 
tion of the vertex correction due to the impurity scatter- 
ing, does not vanish, in general, and, e.g., in 2D k-cubic 
Rashba system [l^ , 2D k-cubic wurtzite system [l^l , and 
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2D k-cubic Dresselhaus system [l^l- Very recently, the 
large SHE in Pt metal at room temperature [l^ was also 
theoretically investigated and was attributed to be an in- 
trinsic one due to the band anti-crossings near the Fermi 
level at the L and X symmetry points in the Brillouin 
zone 

The spin precession around the effective magnetic field 
caused by spin-orbit coupling leads to the fundamental 
problem that the conventionally defined intuitive spin 
current operator i{v, Sz} is not conserved. Therefore, 
how to properly define the spin curren t op erator has been 
intensively studied in recent vear s [2p . l23l .l23.[25l|. In view 
of the spin continuity equation [26l| -I- V • — 
Shi et al. 22] recently provided a proper definition of 
conserved spin current to resolve this issue. The effec- 
tive conserved spin current ^(xsj.) constructed from the 
spin continuity equation is composed of two terms. One 
term is the conventional intuitive spin current operator 
^Sz, and the other term which is so-called torque 
spin current comes from the spin precessional motion. 
Zhang et al.^^ considered the spin Hall coefficients for 
three widely studied semiconductor models, namely, 2D 
k-linear Rashba, 2D k-cubic Rashba and 3D Luttinger 
models, in the clean limit, and found that the conserved 
spin Hall conductivities are dramatically different from 
the conventional spin Hall conductivities. For example, 
in the 2D systems, the conserved spin Hall conductivity 
is equal to the conventional spin Hall conductivity in size 
but has an opposite sign. [2^ In Ref. (23, the results of 
calculations taking into account the conserved spin cur- 
rent as well as impurity scattering effect for 2D k-linear 
Rashba and k-cubic Rashba systems are reported. Re- 
cently, we extended the conserved definition of spin cur- 
rent operator and offered a proper definition of the 0AM 
current operator. ^2811 We also found that in 2D Dres- 
selhaus and Rashba-Dresselhaus systems, the conserved 
spin Hall conductivity is equal to the convention a.1 sp in 
Hall conductivity in size but has an opposite sign.l28j 

Clearly, it is important to consider the new defini- 
tion of spin current [13, [H, ll^, and it is of inter- 
est to know the torque and hence conserved spin Hall 
coefficients in other 2D systems. In the present paper, 
therefore, we study the torque, conventional and con- 
served spin-Hall conductivities in all 2D spin-orbit cou- 
pled systems described by a generic effective Hamilto- 
nian [Eq. ([T])] within the frequency-dependent Kubo lin- 
ear response theory. The generic effective Hamiltonian 
covers all common 2D spin-orbit coupled systems used 
in the literature, such as k-linear Rasha, Dresselhaus, 
Rashba-Dresselhaus, k-cubic Dresselhaus and wurtzite- 
type Hamiltonians (Table I) . We find two interesting uni- 
versal relations among the torque, conventional and to- 
tal conserved spin-Hall conductivities. Furthermore, we 
explore possible connections between conventional spin 
current and orbital motion of carriers and identify the ex- 
istence of a conserved hyper-angular momentum in ro- 
tationally invariant 2D spin-orbit coupled systems. The 
conservation of the hyper-angular momentum 1^ would 



lead to a hyper-selection rule which dictates that the 
up(down)-spin state in the sense of (k x ^)zSz would 
be accompanied by the down(up)-OAM state in these 
systems. 

The present paper is organized as follows. In Sec. II we 
define a generic effective Hamiltonian for 2D spin-orbit 
coupled systems and calculate the time evolution of Pauli 
spin and position operators in the Heisenberg picture. In 
Sec. HI we calculate the conventional and torque spin- 
Hall conductivities by using frequency-dependent Kubo 
formulae and also present universal relations between 
these conductivities. In Sec. IV we report our finding 
that there exists a conserved hyper-angular momentum 
Iz in the systems with the cylindrically symmetric en- 
ergy dispersion. We also demonstrate that the existence 
of Iz leads to the hyper selection rule for the conventional 
spin-Hall current. Our conclusions are given in Sec. V. 
Three appendices to this paper outline our derivation of 
the spin continuity equation, a proof of conservation of 
Iz and a proof of the vanishing of the Iz current in 2D 
rotational invariant systems, respectively. 



II. GENERIC MODEL HAMILTONIAN 

The effective Hamiltonian for spin- 1/2 particles can be 
expressed as a linear combination of Pauli matrices a^^ 
Oy and Oz- In 2D systems, we consider the following 
general effective Hamiltonian, 

IlQ=fl+A{\)G,^B{V)Uy. (1) 

where e° — I?k^ I2m, is the single particle kinetic energy, 
and functions A(k) and -B(k) describe the energy disper- 
sion caused by spin-orbit interaction. In general, A(k) 
and i?(k) can be expressed as A(k) = ri^ki -I- rj^kikj + 

rjf^ihkjki H and S(k) = rifki + rjfjkikj + ijf^ikikjki + 

• • • , where 77's are the coefficients to be determined for 
each specific system. The Einstein summation conven- 
tion is used. The general properties of coefficients ry's are 
determined by the symmetry requirements. For instance, 
time reversal invariance of spin current = "^^^^'^ re- 
quires that A(k) and B(k) must be an odd function of k, 
i.e. A{-k) = -A{k) and B{-k) = -B{k). This leads 
to the fact that the spin dependent part of the Hamil- 
tonian has no spatial inversion symmetry. In appendix 
A, we show that the systems described by Eq. ([T]) sat- 
isfy the spin continuity equation: V • Jj = 7^. In 
Table U we list the specific functions A{k) and B{k) for 
several common 2D systems. However, we should stress 
here that the following derivation is independent of the 
detailed forms of A{k) and B{k). 

We should emphasize that Eq. (1) is an effective 
Hamiltonian for 2D systems valid only near the Bril- 
louin zone center, and is not a bare Hamiltonian that 
describes the band structure of the whole Brillouin zone. 
In other words, Eq. (1) is applicable to the 2D semicon- 
ductor structures with the electron or hole pocket cen- 
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tered at the Brillouin zone center such as p-type zinc- 
blende semiconductors and n-type wurtzite nitrides, but 
not to the metals with a complex Fermi surface such as 
platinum f2l| . In writing the effective Hamiltonian Eq. 
(1), we made the assumption that the particle spin-1/2 
(or the pseudospin-1/2 for k-cubic Rashba Hamiltonian) 
lies in the two dimensional plane. For these spin-1/2 
particles, we need only the two component Bloch wave 
function and thus the effective Hamiltonian can be writ- 
ten as the linear combination of Pauli matrices. Since 
the particle spin lies in the plane, the spin splittings in- 
duced by bulk or structure inversion asymmetry can be 
described by introducing the in-plane components of the 
fc-dependent effective magnetic field, in which they are 
yl(k) and B(k). The periodic potential and spin-orbit 
coupling effect would enter the effective Hamiltonian via 
the A{k) and B{k). The explicit forms of A(k) and B(k) 
depend on the symmetries of the underlying crystalline 
structure and band structure near the Brillouin zone cen- 
ter. Since the 2D system we considered is time reversal 
invariant (zero magnetic field), the spin-splitting would 
result from the spatial inversion asymmetry (or structure 
inversion asymmetry). This implies that A(k) and -B(k) 
are odd functions of k. 

For the convenience of derivation, it turns out to be 
useful to introduce a vector M = {Mx,My). The in- 
plane components of vector M = {Mx,My) are = 
cos 6 = ^ and My = sin 6 — ^. The Hamiltonian [Eq. 
P])] can now be rewritten as 

i/o = efc + A(a X M)„ 



(2) 



where A(k) = {A^ + B^Y^^ is the energy dispersion of 
spin-splitting determined by the explicit forms of A(k) 
and B{k) (e.g.. Table |T|. The vector Pauli matrix used 
in Eq. ^ is a — {(Jx,cry)- The eigenenergy of Eq. ([2]) is 
En{li) = e° — nA(k) and the corresponding eigenvector 
is given by 



-ie(k) 



where the 0(k) is 



m ^ ...... 



(3) 



(4) 



and the band index is denoted as n = ±. It is straightfor- 
ward to show that {axM)l = {Ml + M"^) = 1. The time 
evolution operator exp{iHot/h) can be further written as 



JHot/h 



j4t/h 



nt 



+ i{c7 X M)z sin 



nt 



(5) 

where ft — 2A/h. By using the definition of Heisen- 
berg picture for Schrodinger operator O, 0{t) = 
expiiHot / h)Oexp(—iHQt / h) , one can show that the time 



evolution of Pauli spin operators are given by 
cr^(i) = a^- Mx sin(f7t)(T2 + M^ia ■ M)[cos{nt) ~ 1] 
C7y{t) = (Jx- My s\n{VLt)az + My{a ■ M)[cos{Vlt) - 1] 
cr^(i) = cos{nt)a^_ + {d ■ M) sin(f7i). 

(6) 

It can be shown that (nk|((7 • M)|nk) = (nklcr^lnk) = 
by the use of the eigenstate in Eq. ([3]). This means 
that the expectation value of the z-component of the spin 
operator vanishes in the absence of electric field. The 
time evolution position operator can be written as x(i) = 
x(0) + (5x(t) and 



(5x(t) 



del 



1 



l^+o(^xAf)J — 



hdk ' 2 

1 89 r 

' 2dk 



dk 



t 



{cos{Vlt) - l)a^ + [a ■ M) sin(m) 



(7) 



where x(0) is the initial condition and (5x(t = 0) = 0. 
It can be shown that in the pure Rashba system Eq. 
([7]) would reproduce the result given in Ref. [3(|. The 
physical meaning of each term is as follows. If the spin- 
orbit coupling vanishes, one has A = 0, and the time 
evolution position operator Eq. ^ reduces to the free 

particle equation of motion x(i) = x(0) -I- j-^t- The 
second term of Eq. ([7]) is the displacement arising from 
the anomalous velocity in the presence of spin-orbit cou- 
pling. The anomalous velocity plays an important role 
in the anomalous Hall effect 31]. The third and fourth 
terms have the oscillation behavior inducing the Zitter- 
bewegung [13, [H, [s^ . 



III. SPIN HALL CONDUCTIVITY 

As mentioned before, the conserved spin current is di- 
vided into two terms: 



dt 



(xSz 



dx 
"(ft'' 



^ dt 



(8) 



In addition to the conventional spin current ^Sz, one 
have to introduce the torque spin current x^^ in order 
to satisfy the conserved spin continuity equation. On 
the other hand, the time reversal symmetry of the con- 
served spin current Eq. ([8]) would lead to spatial inver- 
sion asymmetry of spin dependent part of Hamiltonian 
Eq. ([1]). This can be seen as follows. From the com- 
mutator ■^[o'z: Hq\ = f2(T • M, since the position operator 
is even under time reversal operation, the invariance of 
torque spin current under time reversal symmetry must 
require ^(-k) = -A{k) and B{-k) = -J5(k). The time 
evolution of the conserved spin current is 

J{t) - ^{V(t), Sz{t)} + i{x(t), Ho](t)} 

= J^-(t) + J"-(t) + J"«(i), 
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TABLE I: Some common 2D systems where the effective ffamiltonian can be described by Eq. ([T}. The A(k) describes 
the energy dispersion in the presence of spin-orbit coupling, wherein 7((^) = \J cP- fP- — 2al3 sin(2(^), ii{(l>) = i sin(2(;/>), and 
tan~^ (/) — The aJj(O) is the conventionally defined spin-Hall conductivity. The pseudospin angular momentum of the 
k-cubic Rashba hole system used in the calculation of spin current is S = ^ha. Superscript * denotes that the system is not 
rotationally invariant. 



2-D system 


A(k) 


5(k) 


A(k) 




References 


Rashba 


aky 


akx 


ak 


-|e| 

Rtt 


[35] 


Dresselhaus ([001]) 


Pkx 


Pky 


(5k 


e 


[36j 


Dresselhaus ([110])* 


pkx 


-pkx 


\f2pk cos 





[3J 


Rashba-Dresselhaus 




rvA^T — fih,, 


k'yi (n\ 


~\^\ ■ / 2 nl-. 


[28 29 371 

Uin' 


k-cubic Rashba (hole) 


'-^{ki-kX) 


^(kt+kl) 




-9|e|ft2 ^1 1 ^ 
IG-jr'^mait kp kp 


[38J 


k-cubic Dresselhaus* 






/3Dfc^K(9i) 


jel^i"^ f cac4> csc4> ■. 
IG-K^mPo J k+{<f,) k-{(j,y 


[20] 


Wurtzite type 


{OLo + f3ok^)ky 


{ao + pok'^)kx 


aok + Pok? 


tan-i(^™^^^) 
Vom-K V«o/3o 


[10, 39] 



where 



where 



J^^W = ^{v(t), ..(*)} 



is the conventional spin current 



i{5x(t),l[s3,i7o](t) 



(10) 



(11) 



is the torque spin current which is independent of the 
choice of origin of coordinate system and 



J™(t) = i jx(0),l[.s„i/o](i) 



(12) 



is the other part of torque spin current which depends 
on the initial choice of the origin of the coordinate sys- 
tem. The time dependent part of the position opera- 
tor 5x(t) is given by Eq. ([7]). It can be shown that 
{n\.\y-{t = 0)|nk) = (nk|J^-(i = 0)|nk) = and 
(rik|J™(t = 0)|nk) = 0. This leads to the fact that 
the conserved spin current vanishes at t = as required, 
namely, (nk|J'(t = 0)|nk) in the absence of external 
electric field. The value of Eq. is the torque spin 

current with reference to the initial choice of the origin of 
the coordinate system. We could choose the initial posi- 
tion of the carrier as the origin of the coordinate system, 
and as a result, Eq. ([T^ would not contribute to the spin 
accumulation. In that sense, the conserved spin current 
J{t) can be divided into two terms J{t) = J[t) -I- J'^", 



(13) 



corresponds to the total spin current which is free 
from the choice of the origin of the coordinate system. 
Eq. (|13p could satisfy the initial condition, namely, 
(nk|J'(t — 0)|nk) — because it can be shown that 
(nk|J''-(f = 0)|nk) = and (nk|J^-(t = 0)|nk) = 0. 
The frequency-dependent Kubo formula for a spatially 
homogeneous electric field [13] is 



(w i^) Jo 



^E/"k(rik|[J^(t),z;,(0)]|nk) 



(14) 



where q is the carrier charge, i.e., q = — |e| for electrons, 
and /„k is the Fermi distribution at zero temperature. 
The parameter r\ is used to regularize the integral and 
the direction of external electric field is denoted as index 
V. We will calculate the conventional and torque spin- 
Hall conductivities by using Eq. ([H)) . We assume that 
the electric field is applied in the y direction (y — > y). 
The transverse spin current is composed of conventional 
and torque spin-Hall currents, Jx{t^ = Jx'{t) + Jx^i^)- 
The conventional spin-Hall current in the x-direction is 
Jx'i^) = ^{vx,Sz}{t), and it can be evaluated as 



(15) 
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where Vx is defined as Vx = hkx/m. The torque 
spin-Hall current in the x direction is J^' (t) — 
^^{Sx{t), jg[az, Hq]}. After substitution of the commu- 
tator jf-[az,Ho] = Qa ■ M to the torque current and 
straightforward calculation, one can obtain 



do 

{ixt)VLa{t) ■ M + - — sin(m) 



(16) 



where 5x{t) given in Eq. ([7]) was used and $(1) = 
{(Tx{t),ay{ty) wherein (Jx{t) and <Jy{t) are given in Eq. 
([6]). With the definition of conserved spin current, the 
total spin Hall conductivity is the sum of contributions 
of conventional and torque spin-Hall currents, 



(17) 



where Cj = oj + irj. The first and second terms in the 
right hand side of equality correspond to the conventional 
spin Hall current and spin torque current, respectively. 



By using Eq. {H]) and Vy{0) 
-A{a ■ M), one can obtain 



OA 



{a X M), 



hdky 



4 A{kx-^) 



fc — 



(18) 



for the conventional spin-Hall conductivity and 



-2q 



(19) 



for torque spin-Hall conductivity, where Co = lo + irj, 17 = 
2A/?i, is Fermi momentum for band n = ± and dSk — 
kdkd(f>. In the absence of spin-orbit coupling, A = 0, the 
spin Hall conductivity vanishes as one can see from Eq. 
(|T5| and Eq. HH). 



A. The static limit: uj = 

In the intrinsic case and the static limit (i.e. lo — 0), 
after some algebraic calculations, one can obtain from 
Eq. p8|) and Eq. p9|) an exact relation: 



<;(0) = -2<;(0). 



(20) 



Thus, the torque spin-Hall conductivity is simply a con- 
stant (-2) multiple of the conventional spin-Hall con- 
ductivity. This universal relation implies that in the 
spin-1/2 2D systems with spin-orbit coupling the magni- 
tude of the torque spin-Hall conductivity is always twice 
larger than the conventional spin-Hall conductivity and 
has an opposite sign. The spin z-component is not a 
constant of motion, as shown by the commutator of the 
Sz and Hq. We find that, by virtue of the commutation 



properties of Pauli spin-1/2 matrices, [sz,Ho] is a linear 
combination of the in-plane spin components, namely, x 
and y components. We can rewrite the Hamiltonian as 
Hq = + (7 • Beff, where Beff is the effective mag- 
netic field in the k-space. The time derivative of spin 
z-component is daz/dt — (— 2)((t x Beff)z/fi and the 
quantity —2 on the right hand side of equality actually 
yields the result that the magnitude of intrinsic torque 
spin-Hall conductivity is always twice larger than con- 
ventional spin-Hall conductivity. Physically, it is the 
spin precession that leads to the result that the mag- 
nitude of torque spin Hall conductivity is twice larger 
than conventional, and the sign of torque spin Hall con- 
ductivity is opposite to conventional spin Hall conductiv- 
ity. We notice that the k-cubic Rashba system describing 
spin-3/2 heavy hole also obeys Eq. (^0)1 . In the k-cubic 
Rashba system [sl] , Pauli matrices operate on the states 
with spin-3/2 projection along the growth direction. In 
that sense, the k-cubic Rashba system actually repre- 
sents a pseudospin-1/2 system. We also note that the 
spin-dependent part of the Hamiltonian of the k-cubic 
Rashba is originally written as {k^a+ — k'^a-), where 
k± — kx i: iky and (7± — ax i iCy- This can be rewrit- 
ten as A(k)(Tx - B(k)(Ty with A = '-^{k^ - kl) and 
B = ^(fcii + k'!^). The simple relation between intrin- 
sic conventional and torque spin Hall conductivities Eq. 
(PO)) is independent of the detailed forms of spin splitting 
(i.e. A(k) and -B(k)) and henece the energy dispersion 
(i.e. A(k)). The total intrinsic spin Hall conductivity 
^^,(0) = [KyiO] + <l,(0)] is then given by 



'xy 



(0). 



(21) 



The sign of the total spin-Hall conductivity is always op- 
posite to the sign of the conventional spin-Hall conduc- 
tivity. Even if higher order k terms were included in the 
theoretical calculations, the conclusions described above 
would still be true. It should be emphasized that the 
validity of Eq. ((2T|) is independent of the sign of carrier 
charge. Interestingly, Eq. (pTjl suggests that the prop- 
erties of the total intrinsic spin-Hall conductivity can be 
characterized by the conventional spin Hall conductiv- 
ity only. Both Eq. (pO)l and Eq. (|2T|l are the universal 
results for 2D spin-orbit coupled systems. Let us now 
apply formulae Eq. ([TS]) and Eq. to some specific 

systems. In the Rashba-Dresselliaus system, for exam- 
ple, we can derive from Eqs. (fT8)l and (|19p that a^y — 
{q/8T:)sign{a^ — /3^) and ci^^ = —{q/A'K)sign{a'^ — 0^). 
The results agree with previous works for the Rashba sys- 
tem [23^ and for the Rashba-Dresselhaus system ^28j . We 
note that when a — (3, A — —B and 9 is independent of k. 
In this case, the spin-Hall conductivity vanishes [H, [33]- 
For the Dresselhaus Hamiltonian along the [110] direction 
(see Table. I.), we have A = —B and thus the spin-Hall 
conductivity also vanishes. It has been shown that the 
[110] Dresselhaus Hamiltonian with i-z along [001] direc- 
tion is different from the Rashba-Dresselhaus model with 
a = P only by a unitary transformation [40| . The exact 
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SU(2) spin rotation symmetry has been investigated in 
both systems (40|. 

We also note that Eq. pTjl is still true even if there 
exists an uniform in-plane external magnetic field B, as 
can be seen as follows. The interaction of the spin and 
magnetic field is proportional to a ■ B — a^Bx — o-y{— By). 
The Hamiltonian is now given hy H = + <Jx{A{k) + 
I^bBx) — '^y{B{k) — ^isBy), where /^b is the Bohr mag- 
neton. We can redefine functions A'(k) and B'{}s) as 
yl'(k) = A{k) + fiBx and B'(k) = B{k) - fiBy, respec- 
tively. Therefore, though the numerical values of Eq. 
(fTS]) and Eq. HH) may change, Eq. ^ is still valid in 
the presence of an uniform in-plane magnetic field. 

Our predictions of, e.g., the conserved spin Hall cur- 
rent and conductivity [Eq. (21)], can be tested by direct 
measurements of the spin Hall current or conductivity. 
In particular, our prediction that the total spin-Hall con- 
ductivity differs from the conventional spin-Hall conduc- 
tivity only in sign [Eq. (21)], could be easily tested. As 
for the induction of a magnetic field by a charge current, 
a spin-current would generate an electric field [gj, [4l| . 
Therefore, our prediction for the conserved spin Hall cur- 
rent could be tested by comparison of the measured spin 
current-induced electric field to the numerical simulation 
for, e.g., a mesoscopic spin-orbit coupled system. An- 
other kind of experiments is to determine the inverse spin 
Hall conductivity by measuring the charge current and 
Zeeman field gradient [12, [11] . The sign and magnitude 
of the conserved spin Hall conductivity can then be ob- 
tained via the Onsager relation [l^, [1^. Indeed, large 
spin Hall effect in metallic systems at room temperature 
has recentl y b een detected by the method of inverse spin 
Hall effect hj, [l^ ITsj . Hopefully, our interesting predic- 
tions would stimulate measurements of inverse spin-Hall 
effect in 2D semiconductor systems in the near future. 

B. Finite frequency case: Cj 

At finite frequencies, the torque spin-Hall conductivity 
is not a constant multiple of conventional spin-Hall con- 
ductivity. Nonetheless, we find that Eq. pS]) is related 
to Eq. ([TO)) by the following equation: 

alliCo)^i-2-u-^)a^/y{u). (22) 

The second term in Eq. (j22p comes from the variation 
of the k-space effective magnetic field with the frequency 
dependent external electric field. Taking into account the 
torque spin current, we find that the total spin-Hall con- 
ductivity is directly related to the conventional spin-Hall 
conductivity still. Substituting Eq. into Eq. (fTT]) . 

we obtain the relationship between conventional spin- 
Hall conductivity fT^|^(a)) and total spin Hall conductivity 
axy{uj) in the presence of non-zero frequency-dependent 
electric field: 

= l^'^ili^)] ■ (23) 



Eq. ([23)1 shows that the total spin-Hall conductivity can 
be determined directly from the frequency spectrum of 
conventional spin-Hall conductivity. Unlike the static 
limit, the total ac spin-Hall conductivity is not pro- 
portional to the conventional spin-Hall conductivity. It 
follows from Eq. ([23|) that the step function behavior 
of the conventional spin-Hall conductivity would result 
in a large response of the total spin-Hall conductivity. 
This large response has recently been investigated in the 
Rashba-Dresselhaus system in Ref. [l^ It can be shown 
that Eqs. (HH]) and (fTOj) for the Rashba-Dresselhaus sys- 
tem agree with the results in Ref. [2^. We also find that 
the simple relation between the torque and conventional 
spin-Hall conductivities [Eq. ([22]) ] would be maintained 
even when the external magnetic field is applied. 

IV. CONSTANT OF MOTION IN 
ROTATIONALLY INVARIANT SYSTEMS 

A. Constant of motion 

If a 2D system is invariant under rotation about the 
z-axis, its energy dispersion is cylindrically symmetric, 
i.e., A = J2e (^ = 1, 2, 3, • • • ). We find that in this 
case, there exists a conserved quantity whose operator 
is defined as 

= (kx 1^) ^z+Lz- (24) 

Let us call this quantity Iz the hyper-angular momentum. 
We can show that Iz satisfies the following commutation 
relation (see appendix B): 

[Iz,Ho] = 0, (25) 

where Sz — f fz, i-^z = fi(xky — yk^) is the z-component 
of the 0AM, and Hq is given in Eq. (P) or Eq. 
Interestingly, this implies that in the rotationally invari- 
ant spin-orbit coupled systems, the flow of (k x f^)^ Sz 
would be accompanied by the orbital angular momentum 
Lz, and the combination of these quantities is actually a 
constant of motion. As a result, the spin current would in 
general be accompanied by the 0AM current because of 
the spin-orbit coupling. In that sense, the external elec- 
tric field would induce the current of angular momentum 
(k X )zSz and the 0AM current simultaneously. 

Further, it can be shown that the hyper-angular mo- 
mentum Hall current {^{Iz,Vx}) vanishes in the steady 
state case within the linear-response Kubo formalism (see 
appendix C). It follows that the up-spin (down-spin) 
state in the sense of (k x f^)zSz would be accompa- 
nied by the down-OAM (up-OAM) state, rendering the 
hyper-angular momentum conserved. Let us now apply 
this result to some specific systems. For the wurtzite 
type and k-linear Rashba systems, Iz — Sz + Lz, i-e., the 
hyper-angular momentum is equal to the total angular 
momentum. For the k-cubic Rashba hole system, one 
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would have Iz = 3s z + Lz- Therefore, Eq. (p4)) gives 
the correct pseudospin angular momentum of hole which 
is |S(Tz. The conservation of hyper- angular momentum 
in these systems would then lead to the result that the 
spatial trajectory of a down-spin {—iz) carrier would be- 
have as having its orbital angular momentum pointed to 
+ez and vice versa (see Fig. 1. (a)). In other words, 
in these systems, the situation of an up-spin (down-spin) 
state accompanied by an up-OAM (down-OAM) trajec- 
tory is forbidden. This is a hyper-selection rule that is 
present in the cylindrically symmetric 2D spin-orbit cou- 
pled systems. It must be emphasized that this hyper 
selection rule depends on the quantity (k x |f )z. For ex- 
ample, in the k-linear Dresselhaus system, we find that 
Iz — —Sz + Lz- Therefore, in contrast to the k-linear 
and wurtzite type systems, the hyper-selection rule in 
the k-linear Dresselhaus system implies that an up-spin 
(down-spin) state would be accompanied by an up-OAM 
(down-OAM) state [see Fig. 1. (b)]. The situation of 
an up-spin state accompanied by an down-OAM is now 
forbidden in this Dresselhaus system. Nevertheless, we 
should emphasize that the electric field-induced 0AM 
current would not result in magnetization accumulation 
at the edges of sample. This is due to the fact that the 
0AM is not an intrinsic quantity of electrons or holes, 
i.e., the magnetic moment associated with the orbital an- 
gular momentum would vanish when the carrier velocity 
reaches zero at the edges of sample, as can be understood 
from the definition of Lz- Therefore, the magnetic mo- 
ment accumulation at the edges of sample would come 
from the spin angular momentum only. In short, it is in- 
teresting to notice that the topological quantity {k.x ^)z 
is a integer number for k-linear Rashba, k-linear Dressel- 
haus, k-cubic Rashba and wurtzite-type system, whereas 
they are 1, -1, 3, 1, respectively. The topological number 
indicates that the constant of motion in 2D rotationally 
invariant system is the hyper-angular momentum rather 
than simply the total angular momentum Sz + Lz- 

It should be pointed out that although both hyper- 
angular momentum conservation and hyper-selection rule 
exist in rotationally invariant systems, this rotational 
symmetry may be broken when higher order terms of 
ki are included in the A(k) and B(k) in Eq. (1). It 
would be necessary to include the higher order terms of 
ki in Eq. (1) when the 2D semiconductor systems consid- 
ered have a very large carrier concentration. When the 
weak symmetry-breaking higher order terms do appear, 
the hyper-angular momentum [Eq. (24)] is no longer con- 
served but the universal relation [Eq. (20)] still holds. 

In the next subsection, we will describe the close re- 
lation between Berry vector potential and hyper-angular 
momentum. 



B. Geometrical interpretation 

As mentioned above, the quantity (kx ^)zSz together 
with the orbital angular momentum is conserved in a ro- 



a) 



up-OAM 



up-spin 



A- 



down-spin 



down-OAM 



b) 



up-OAM 



down-spin 



r" 



down-OAM 



up-spin 



FIG. 1: Schematic diagrams showing the relationship between 
spin and orbital angular momentum (0AM) of a carrier. The 
direction of the OAM is determined by the right-handed sense. 
The dotted line illustrates the spatial trajectory of a carrier, 
(a) for the Wurtzite type, k-linear and k-cubic Rashba systems 
and (b) for the k-linear Dresselhaus system. 




FIG. 2: A schematic diagram showing the relationship be- 



tween vectors k. 



ak' 



and k x 



tationally invariant system. In the following, let us ex- 
plain that the quantity (k x ^)z actually comes from 
the topological properties of the Berry vector potential. 
In the 2D systems with a cylindrically symmetric disper- 
sion, the vector in general, is perpendicular to the 
wave vector k, i.e., k • ^ =0. For example, it can be 
shown that the wurtzite type, k-linear Rashba, k-cubic 
Rashba and k-linear Dresselhaus systems, the dot prod- 
uct of k and ^ is zero. In that sense, the three vectors 
k, ^ and k x form an orthogonal frame fixed on the 
carrier (see Fig. 2.). 

On the other hand, the Berry vector potential in a 
system described by Hamiltonian Eq. ^ can be written 
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^(k) 



{nk\i—\nk) 



1 de 

2dk' 



(26) 



where the eigenvector Eq. ([3]) was used. The projec- 
tion of spin operator onto the in-plane axes contains two 
terms. One term is the so-called spin helicity ct • k in 
the 2D spin-orbit coupled system, and the other is the 
projection of spin on the Berry vector potential ■ a. 
The quantity (k x ^)zSz then comes from the noncom- 
mutativeness of the two in-plane projections, viz. 



(27) 



where the commutation relations of Pauli matrices were 
used. The overall coefficient of (k x ^)zSz cannot be 
determined by the commutation relation alone. How- 
ever, in the system with the cylindrically symmetric 
dispersion, the hyper-angular momentum conservation 
forces the overall coefficient of (k x ^)zS2 to be unity. 
Similar to the topological force induced by the non- 
commutative position operator [4], the spin part of con- 
servation of the hyper-angular momentum comes from 
the non-commutative properties oi a ■ .s^f and a ■ k, 
whereas they are the projection of spin on the two or- 
thogonal axes. Finally, it must be stressed that in the 
systems with the non-cylindrically symmetric dispersion, 
the quantity (k x ^)zSz can also be defined as the non- 
commutativeness of spin helicity (j-k and a-si . However, 
in that case, the quantity (kx ^^)zSz plus the orbital an- 
gular momentum is not conserved. Therefore, the orbital 
motion of carrier does not accompany with (k x )zS2, 
namely, the hyper-angular angular momentum is not con- 
served in this case. 

The conserved quantity in spin-orbit coupled systems 
can be written as the sum of the spin and orbital terms. 
In the free atomic case, the spin term is just the Pauli 
spin operator. When the crystal environment is included, 
the spin term appears to be different from the Pauli spin 
operator. It has a non-trivial dependence on energy dis- 
persion that arises from the spin-orbit coupled effect. We 
find that the general coefficient is a non-trivial multi- 
plication of (k X f^)z in rotationally invariant system. 
The quantity (k x ^)z is proportional to the expecta- 
tion value of the orbital angular momentum. 

It is of course not surprising that a conserved angular 
momentum would exist in rotationally invariant systems. 
The conserved angular momentum in spin-orbit coupled 
systems would be the sum of the spin and orbital terms. 
In the free atom case, the spin term is just the Pauli spin 
operator, and the sum of the spin and orbital terms is in- 
deed the total angular momentum. However, we find here 
that in the presence of the crystal environment, the spin 
term is not necessarily equal to the Pauli spin operator, 
but has a non-trivial dependence on the energy disper- 
sion [Eq. (24)] that arises from the spin-orbit coupling 



effect. Therefore, we would like to use the hyper-angular 
momentum here to differentiate Eq. (24) from the well- 
known expression of the total angular momentum in the 
free atom case. Furthermore, in the systems described by 
the generic Hamiltonian Eq. (1), we find that the coeffi- 
cient in the spin term in Eq. (24) is related to the Berry 
vector potential of the underlying band structure. 



V. CONCLUSIONS 

In conclusion, we have derived some interesting re- 
lations between the conventional and torque spin-Hall 
conductivities for all 2D spin-orbit coupled systems de- 
scribed by the generic effective Hamiltonian Eq. ([1]) in 
the presence of frequency-dependent electric field. In 
particular, we find that an universal relation [Eq. (I2ip ] 
between total and conventional spin-Hall conductivities, 
i.e., cr|y(0) = — cr^^(O) for w = 0. Eq. (PT|) is independent 
of the detailed form of energy dispersion A(k) (i.e. A(k) 
and S(k)), and hence its validity is not restricted to the 
systems listed in Table. HI We also found that in 2D ro- 
tationally invariant systems, a conserved hyper-angular 
momentum exists and the hyper-angular momentum 
current vanishes. This would result in a hyper-selection 
rule that the up-spin (down-spin) state in the sense of 
(k X '^)zSz would be accompanied by the down-OAM 
(up-OAM) state. Finally, we explained that the spin 
dependent part of Iz comes from the noncommutative 
property of spin helicity u ■ k and a ■ . 
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APPENDIX A: SPIN CONTINUITY EQUATION 

In this appendix, we show that the Hamiltonian 

H=^+E,, + F(x) (Al) 
2m 

would satisfy the spin continuity equation: 

^+V-3^=Tz, (A2) 

where Hgo = ^(p)^-,; — B{p)ay describes the spin-orbit 
interaction, V{x) = — gE • x is the potential induced 
by a homogeneous electric field E, Sz — '^^Sz'^ is the 
spin density, J* — Re ["ii^ ^{v, Sz}'^] is the conventional 
spin current and 7^ — Re jjr[sz, Hq]'^] is the source 
term of spin current. The carrier velocity is defined as 
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V = The real part (imaginary part) of [• • •] is de- 
noted as Re[- ■ ■] {Im[---]), and Sz — ^(Tz- By using 
Schrodinger equation with the two-component wave func- 



tion = 



one can obtain: 
d 

ih—i'^hz'i') 



0* 

ih— = H^, (A3) 



-{^M,)Ksz^ + ^^Szi^^) 

Zm Zm 

+ {-(y(x)^')ts;,* + *^s2y(x)*} . 

(A4) 



We note that [(i/^o + F(x))^']ts^* = [*ts^(H^^ + 
y(x))^']'l' because the Pauh spin matrix az satisfies o-| = 
(Tzii.e., is Hermitian. On the other hand, one can also 
show that 

[^^s^p^* - (p2'J')^s^*] = 2p • i?e[^'tps^*]. (A5) 
After substituting Eq. (jA5P into Eq. (jA4|, one gets 
d 



ft— (^-ts^*) =p • i?e[*tZs^^] 



dt 

+ 2iIm[^i^SzHso'i'] + 2iIm[^i^SzV{-x.)'^]. 

(A6) 

We also note that SzHso can be written as SzHso = 
^{Sz, 7?so} -|- ^[sz, ^^so], where the first term vanishes be- 
cause the Pauli matrices satisfy {(7^, cTj} = 25ij. One can 
obtain 



ft 



-Re 



^^-[Sz^H.o]'^ 

in 



(A7) 



The last term of Eq. (jA6p vanishes because T^(x) is real 
as required by the Hermitian property of Hamiltonian 
Eq. Finally, substitution of Eq. ((XtI) into Eq. 

dMI yields 



-V -Re 
+ Re 



«ft 



(A8) 



where the commutation relation { gp" i ■^z} = 
used. Eq. (|A8|) is the desired spin continuity equa- 
tion. The average spin torque vanishes [11], and hence 
we have / dxT^ = 0. The spin torque density can be 
written as the divergence of spin torque dipole density 
Pr(x), namely, 7^ = — V • Pr(x). On the other hand, 
the spin dipole density vanishes outside the sample, and 



we have dxP^ = dx(— xV ■ Pt) = /ydx(x7^). 
Therefore, the spin dipole density can be written as 
P^(x) = Re[^myi,^}^]. Finally, the effective con- 
served spin continuity equation can be written as 



dSz 
dt 



(A9) 



where J'(x) = i?e[5'^J'^'] and the effective conserved 
spin current operator is J — ^{v, Sz} + |{x, ^} which 
is the sum of conventional and torque spin currents. 



APPENDIX B: CONSERVATION OF h 

In this appendix, we demonstrate that the hyper- 
angular momentum Iz defined in Eq. (|24p is a conserved 
quantity when the energy dispersion A is rotationally in- 
variant. First, it can be shown that the velocity operator 
can be written as 



1 

ih 



[X, Hq] = Wx + J-g^i'^ X M)z + 



09 

hdk„ 



•A(ct • M), 
(Bl) 

where Hq = e° -t- A((7 x AI)z was used. The y-component 
of velocity can be obtained by replacing the index x by 
y. The z-component of orbital angular momentum oper- 
ator is defined as Lz = h{xky — ykx)- Using the velocity 
operator, the commutator [Lz^Hq] is straightforwardly 
evaluated as follows: 



86 



00 



aA, aA, 



(B2) 



We now define the operator Iz as 



Sz + Lz 



(B3) 



where 9 ~ tan ^ (5) ^ind Sz — ^ ctz is the Pauli spin 
operator. Using Eqs. (|B2[) and (jBSp . we obtain 



1 /9A aA 



ih 



Oky 



{ax M)z, 



(B4) 



where -^[az^Ho] = ■ M was used. Eq. (|B4)l is the 
main result of this appendix, and it means that, in gen- 
eral, the Iz operator is not a conserved quantity. The 
right hand side of Eq. (jB4p explicitly depends on the 
form of energy dispersion. It is interesting to note that 
the spin term of hyper- angular momentum (Eq. (|B3|) ) 
is not the Pauli matrices with the multiplication of ft/2, 
but the multiplication of 



1 



(B5) 
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which is further explained in Sec. IV. B. 

Now consider a 2D system with the cyhndrically sym- 
metric energy dispersion that can be written as the power 
series of magnitude of k denoted as k = |k|, namely, 
A = ^^Cf/c^. The right hand side of equahty in Eq. 
()B4p then becomes 

(B6) 

Therefore, the hyper-angular momentum is a con- 
served quantity in the rotationally invariant 2D systems. 



APPENDIX C: NULL HYPER- ANGULAR 
MOMENTUM CURRENT 

In this appendix, we will show that the hyper-angular 
momentum current ^{Iz,'v} vanishes in the linear re- 
sponse regime. In the static case, the Kubo formula can 
be written as 



/nk fn'li. 



(i?„(k)-ii;„,(k))2 



X Im(nk| J^|n'k)(n'k|z;t,|nk) 



(Cl) 



where both eigenstate \nk) and eigenenergy i5„(k) are 
given in Sec II. In the following, the external electric 
field is assumed to be applied in the y direction, and we 
calculate the conductivity axy First of all, the hyper- 
angular momentum current can be divided into two terms 
J^^ = J'^' + J^' . First term is the hyper-spin current 



J: 



^{(kx^^) (C2) 



corresponding to the hyper spin-Hall conductivity ct;^^ 
The second term of jj^ is the orbital current 



(C3) 



corresponding to the orbital-Hall conductivity <J^y- 

a result, the hyper angular momentum Hall conductivity 

can be written as 



xy xy ' xy " 



(C4) 



We first calculate the orbital-Hall conductivity. By using 
the velocity operator Eq. (jBip . we have 



(C5) 



Im{nk\J^' \ - nk){-nk\vy\nk) 

hkx ^ de fde ■ 

X rC 



^ 2m dky V 9k 



de ,.dA dB. _^g^ fde , 

)e } I — X k 



Yn'^^^'^dk^'dk^^dk.. 



\dk 



The second term of Eq. (|C5|1 does not contribute to the 
orbital-Hall conductivity a^^ because of the even power 
of band index n. The orbital-Hall conductivity cr^^ with 



substitution of Eq. (IC5|) gives 
qh 



2 /■'^^ fc-^ (de 



(C6) 



We now consider the hyper-spin Hall current. Taking 
into account the hyper-spin and the velocity operator Eq. 
(|B1[) . we have 



Im{nk\jf'\ - nk) {-nk\vy\nk) 

hk, ^ de / de 

= n Att^ I k X 



2m dk,, 



dk 



(C7) 



Inserting Eq. (|C7p into a'^y , we obtain 



qh 



2 k^m: (. de 



Comparison of Eq. (|C6|) and Eq. (|C8|) gives cr^^ = 
—fJ^y- As a result, we have a^^y = a^^ + a^^ = 0, i.e., 
the hyper-angular momentum current is zero. 
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